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基于 高 斯 扰动 和 指数 递减 策略 的 改进 蝙 量 算法 


= 昌 吕 
(1. 郑州 财经 学 院 现代 物流 与 管理 系 ,郑州 450000; 2. 河南 大 学 商学 院 , 河南 开封 475004; 3. 河南 大 学 管理 科 
学 与 工程 研究 所 , 河南 开封 475004) 


摘 要 : 针对 基本 蝙蝠 算 法 后 期 收 钙 速度 慢 、 收 化 精度 不 高 、 稳 定性 不 强 等 问题 ， 提 出 基于 高 斯 扰动 和 指数 递减 策 
略 的 改进 蝙蝠 算法 (GDEDBA)。 把 指数 递减 策略 引入 速度 更 新 公式 ， 使 算法 迅速 进入 局 部 寻 优 并 展开 精确 搜索 ; 
构造 高 斯 扰动 项 加 入 到 局 部 新 解 产 生 公式 ， 使 局 部 新 解 中 所 有 粒子 与 当前 全 局 最 优 粒子 产生 信息 交流 与 学 习 ， 防 止 
陷入 局 部 最 优 ， 增 加 种 群 多 样 性 ; 设计 扰动 控制 因子 来 控制 高 斯 扰动 的 扰动 范围 ， 增 强 算法 的 稳定 性 。15 个 测试 函 
数 的 仿真 结果 表明 ， 改 进 算法 的 寻 优 性 能 显著 提高 ， 收 你 速度 更 快 ， 求 解 精 度 更 高 ， 稳 定性 更 强 。 

关键 词 : 蝙蝠 算法 ; 高 斯 扰动 ; 指数 递减 策略 ; 算法 改进 ; 函数 优化 
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Improved bat algorithm based on Gaussian disturbance and exponential decreasing strategy 


Song Yiminl 2 LiYu2 六 
pe (1. College of Logistics & Management, Zhengzhou Institute of Finance & Economics, Zhengzhou 450000, China; 2. 
= Business School, Henan University, Kaifeng Henan 475004, China; 3. Research Institute of Management Science & 
Engineering, Henan University, Kaifeng Henan 475004, China) 


Abstract: Aiming at the shortcomings of the basic bat algorithm such as slow convergence, low convergence precision and 
weak stability, this paper designed an improved bat algorithm (GDEDBA) based on Gaussian disturbance and exponential 
decreasing strategy. It introduced the exponential decreasing strategy into the speed update formula, could enable the 
algorithm to enter local optimization quickly and exactly; it added the constructed Gaussian disturbance term to the local 
new solution generation formula, then the information exchange and study between all particles in the local new solution 
and the current global optimal particles cloud happened, prevented falling into local optimum and increased population 
diversity; it designed the disturbance control factor to control the disturbance range of Gaussian disturbance, enhanced the 
stability of the algorithm. The simulation results of 15 classical test functions shows that the optimization performance of 
improved algorithm is significantly improved, the convergence speed is faster, the solution accuracy is higher, and the 
stability is stronger. 

Key words: bat algorithm; Gaussian disturbance; exponential decreasing; algorithm improvement; function optimization 


0 引言 一 种 会 飞 的 兽 类 。 其 中 的 微型 蝙蝠 由 于 视觉 退化 ， 靠 回声 定 
二 位 来 代替 视觉 功能 。 微 型 蝙蝠 在 飞行 过 程 中 ， 不 断 发 出 25 一 
2010 年 ， 剑 桥 大 学 学 者 Yang 通过 模拟 蝙蝠 回声 定位 行 ”150 千 赫 效 范 围 内 的 超声 波 ， 超 声波 遇 到 障碍 物 返 回 到 微型 
为 ， 提 出 一 种 新 的 智能 优化 算法 一 一 蝙蝠 算法 (bat-inspired ”蝙蝠 的 双 耳 ， 微 型 蝙蝠 通过 判断 发 出 声波 到 收 到 返回 声波 的 
algorithm)D。 蝙 晤 算法 自 提 出 以 来 ， 以 其 模型 简单 、 易 于 实 。 时 间 差 感知 障碍 物 的 距离 ， 通 过 返回 声波 的 时 间 间 距 感 知 障 
现 等 优点 被 广泛 应 用 于 多 个 领域 ,如 函数 优化 中、 工程 结构 。 碍 物 在 哪里 ， 利 用 回声 的 响 度 变化 感知 障碍 物 是 什么 ， 以 此 
Ob ee A ie ls en ag 构建 三 维 立 体 图 像 。 微 型 蝙蝠 利用 多 普 勒 效应 对 周围 场景 作 
。 但 是 蝙蝠 算法 存在 后 期 收敛 速度 慢 、 易 陷入 局 部 最 优 ， 出 判断 ， 离 猎物 越 近 ， 波 长 越 短 ， 频 率 越 高 ， 离 猎物 越 远 ， 
Dot We de ed oad a 波长 越 大 ， 频 率 越 低 ， 比 并 且 可 以 判断 猎物 的 移动 速度 。 
等 人 [3 用 惯性 权重 因子 来 控制 上 一 步 速度 的 影响 大 小 ， 从 而 蝙蝠 算法 模拟 微型 蝙蝠 搜索 猎物 ， 是 一 种 基于 种 群 的 随 
来 平衡 全 局 和 局 部 搜索 之 间 的 强度 。 Meng 等 人 中 融合 栖息 。 机 寻 优 算法 ， 问 题 的 解 被 抽象 成 搜索 空间 中 的 蝙蝠 粒子 ， 每 
地 选择 和 多 普 勒 效应 改进 蝙蝠 算法 (NBA)， 进 一 步 模 仿 了 蝙 。 只 蝙蝠 都 有 一 个 由 优化 问题 决定 的 适应 度 与 之 对 应 ， 蝙 晤 个 
蝠 的 生物 学 行为 对 蝙蝠 算法 改进 。Luo 等 人 04 提 出 双子 群 莱 栖 通过 调整 频率 、 响 度 、 脉 冲 发 射 率 ， 追 随 当前 的 最 优 蝙 晤 
维 飞行 蝙蝠 算法 (DLBA)。 此 外 还 有 混合 蝙蝠 算法 5~1]。 针 在 解 空间 中 搜索 。 蝙 蝠 粒子 在 位 置 xX; 以 速度 V 随机 飞行 ， 以 
对 蝙 蝙 算 法 的 不 足 , 本 文 提 出 一 种 新 的 改进 算法 (GDEDBA)。 定 频 率 名 (或 14)、 可 变化 波长 14 (或 1) 和 响 度 4 搜索 
和 猎物 ， 根 据 猎 物 与 自己 的 距离 调节 发 射出 的 脉冲 波长 (或 频 
1 。 基本 蝙 量 算 法 率 )， 并 在 靠近 狂 物 时 调整 发 射 脉冲 的 频 度 rs[0U ， 其 中 0 
蝙蝠 是 一 种 神奇 的 物种 ， 是 地 球 上 目前 为 止 发 现 的 唯一 表示 没有 脉冲 ， 表 示 刚 发 现 猫 物 ， 在 接近 猎物 的 过 程 中 ， 逐 
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渐 调整 为 最 大 脉冲 1。 而 响 度 则 在 搜索 猎物 过 程 中 从 最 大 值 
( 正 值 〉 4 到 最 小 值 A。 。 
蝙蝠 位 置 X! 和 速度 Vi’ 在 步骤 1 时 的 更 新 公式 如 下 : 


fi= frin + (fn ~ fain)p ， (1) 
Vi =Vm + XP xX" )f, O) 
X= Xr + G) 


其 中 :ps[0 是 一 个 随机 向 量 ，X 是 当前 全 局 最 优 位 置 。 
微型 蝙蝠 具有 群体 学 习 能 力 ， 能 追随 当前 最 接近 猪 物 的 
蝙蝠 个 体 搜 索 猫 物 。 受 此 启发 ， 蝙 蝠 算法 设计 局 部 新 解 的 更 
新 机 制 ， 在 当前 全 局 最 优 解 的 附近 产生 局 部 新 解 ， 局 部 新 解 
的 更 新 公式 如 下 : 


Xio = X™ teA', ee[-h1] (4) 
其 中 : se[-1 是 随机 数 ;， 4 =< 入 > 是 所 有 蝙蝠 在 这 一 代 里 的 
平均 响 度 。 


脉冲 发 射 的 响 度 4 和 速率 “的 更 新 公式 如 下 : 

汪 =a4 (5) 
= ee (6) 
其 中 : wx 和 7 是 恒 量 ,0< wc <1 和 7 >0。 

2 ”基于 高 斯 扰动 和 指数 递减 策略 的 改进 蝙蝠 算法 
2.1 指数 递减 策略 分 析 

蝙蝠 算法 和 粒子 群 算法 都 是 通过 粒子 在 解 空 间 中 进行 全 
局 搜索 和 局 部 精确 搜索 相 结 合 的 方法 来 实现 最 终 寻 优 的 。 有 
效 的 全 局 搜索 可 以 帮助 粒子 跳出 局 部 最 优 解 ， 有 效 的 局 部 搜 
索 则 帮助 粒子 在 当前 搜索 区 域 进行 精确 搜索 ,利于 算法 收敛 。 
那么 ， 如 何 来 平衡 全 局 搜索 能 力 和 局 部 搜索 能 力 ， 无 疑 会 对 
算法 的 寻 优 性 能 产生 巨大 影响 。 为 此 ， 许 多 专家 学 者 进行 了 
大 量 研究 工作 ，PSO 算法 在 此 方面 的 改进 策略 主要 有 线性 递 
减 权 重 策略 091、 自 适应 权重 策略 Po 、 随 机 权重 策略 P 。 

文献 [22] 受 线性 递减 权重 策略 的 启发 ， 构 造 三 种 不 同 策 
略 对 粒子 群 算法 改进 。 实 验 结果 表明 ， 指 数 递减 策略 效果 最 
佳 。 针 对 基本 蝙蝠 算法 后 期 收敛 速度 慢 且 易 在 全 局 最 优 解 附 
近 产 生 震 荡 的 不 足 ， 受 文献 [22] 启 发 ,用 指数 递减 策略 对 蝙蝠 
算法 进行 改进 。 指 数 递 减 策略 因子 w 的 表达 式 为 
W = Waia (Was / Win ) ro) (7) 
其 中 :wow 、 wm 分 别 代表 w 的 最 大 值 和 最 小 值 ，1 代表 当前 
迭代 次 数 ，tw 代表 最 大 连 代 次 数 ; 参数 < 的 取 值 选取 与 文献 
[22] 一 致 ，c=10， 取 ww =0.9， wmn =0.4。 


一 ; 


0.4 一 开 入 一 
0 100 200 300 400 500 600 700 800 900 1000 
Generations 


图 1 指数 递减 策略 因子 w 的 迭代 曲线 
Fig.1 Iterative curve of exponential decreasing strategy factor w 
指数 递减 策略 因子 随 迭 代 次 数 的 变化 如 图 1 所 示 。 通 过 
观察 : w 在 迭代 初期 随 着 迭代 次 数 增加 以 很 快 的 速率 递减 ， 
而 在 迭代 的 中 后 期 ，w 下降 速度 逐渐 变 得 缓慢 。 在 基本 蝙蝠 
算法 中 ，w 是 个 恒定 值 1， 容 易 造成 算法 在 迭代 前 期 陷入 局 
部 最 优 ， 在 迭代 后 期 收敛 速度 变 慢 。 这 里 将 指数 递减 策略 因 


脉冲 发 射 的 响 度 和 和 速率 也 要 随 着 迭代 过 程 进行 更 新 。 
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子 w 加 入 到 基本 蝙蝠 算法 的 速度 更 新 公式 中 ， 可 以 使 算法 迅 
速 进入 局 部 寻 优 ， 并 有 效 展开 精确 搜索 ， 确 保 找 到 最 优 解 。 
这 样 一 来 ， 既 有 利于 提高 算法 的 收敛 速度 ， 又 有 利于 提高 算 
法 的 寻 优 精度 。 
为 了 测试 该 策略 对 算法 性 能 的 有 效 性 ， 对 比分 析 本 文 的 
改进 算法 (GDEDBA) 与 删除 该 策略 的 算法 (BA1)， 在 保 
证 相关 参数 一 致 的 基础 上 ， 分 别 独立 运行 五 个 测试 函数 ， 取 
其 平均 值 进 行 求解 精度 对 比 ,结果 如 表 1 所 示 。 给 出 GDEDBA 
与 BA1 运行 函数 Fl 的 收敛 对 比 图 进行 收敛 速度 对 比 ， 结 果 

如 图 2 所 示 。 
表 1 GDEDBA 与 BA1 的 仿真 结果 


Table 1 Simulation results of GDEDBA &BA! 
Fl F2 F3 F4 FS 
BAI1l 5.79E-13 2.73E-06 4.42E-14 5.56E-15 0 
GDEDBA 9.87E-121 4.27E-61 0 0 0 
1020 a 
BA1 
i ~ GDEDBA 
jo 
@ 10" 
oo 
后 
号 10™ 
10™ 
10100 
10120 上 


300 400 500 600 700 800 900 1000 
Generations 


0 100 200 


图 2 GDEDBA 与 BA1 运行 Fl 的 收敛 曲线 
Fig.2 Iterative curve of flabout GDEDBA &BAIl 
分 析 表 1 和 2 可 得 ， 删 除 指数 递减 惯性 权重 后 ， 算 法 的 
收 和 敛 速度 和 收敛 精度 大 幅 下 降 ， 说 明 指 数 递减 惯性 权重 可 以 
有 效 提高 算法 的 寻 优 性 能 。 那 么 ， 改 进 后 算法 的 速度 更 新 公 
式 为 


中 三 WwWV +(Xi—X")f (8) 

2.2 高 斯 扰动 过 程 分 析 
观察 基本 蝙蝠 算法 中 局 部 新 解 的 更 新 公式 ( 式 (4))， 可 以 
发 现 它 是 由 当前 全 局 最 优 解 和 随机 数 与 响 度 的 乘积 项 简单 相 
加 构成 ， 新 解 粒子 与 当前 全 局 最 优 粒子 之 间 没 有 产生 信息 交 
流 和 学 习 , 全 靠 响 度 项 来 调节 ,这 就 使 得 对 当前 全 局 最 优 粒子 
的 有 效 信息 利用 能 力 不 高 ， 使 新 解 粒 子 丧 失学 习 能 力 ， 还 会 
带 来 种 群 多 样 性 的 损失 ,也 不 符合 生物 习性 。 
对 此 ， 这 里 在 产生 局 部 新 解 的 公式 Xi。, =X “+<s4 中 ， 对 
前 全 局 最 优 位 置 的 每 一 个 粒子 都 加 入 一 个 高 斯 扰动 项 3 


Ik 


Gaussiantp,o?) ,其 中 : 表示 均值 ; 表示 方差 ; 4=0,07 =|X%|，; 


lx 表示 第 1 次 迭代 时 种 群 的 全 局 最 优 粒 子 的 绝对 值 。 改 进 
后 的 局 部 新 解 产 生 公 式 描述 如 下 : 


Xi = Xta DGaussian(k,o)+eA’ , £ e[-1,1] (9) 


其 中 ;p=0,0? =|X"%|; ss 是 随机 数 ，4 =<4> 是 所 有 蝙 
蝠 在 这 一 代 里 的 平均 响 度 ; a 表示 扰动 控 秆 
对 点 乘法 。 
如 此 ,新 解 的 产生 是 由 当前 全 局 最 优 解 、 均 值 为 0， 方 差 


A 


妹子 ; ，@ 表 示 点 
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宋 一 民 ， 


为 当前 全 局 最 优 粒 子 绝对 值 的 高 斯 扰动 项 以 及 响 度 项 三 项 构 


成 ， 更 加 符合 蝙蝠 之 


点 的 实际 。 


2.3 ”扰动 控制 因子 


可 以 观察 式 (4)， 


基本 蝙蝠 算法 存在 求 


间 会 有 信息 交流 和 相互 学 习 的 


生物 学 特 


识 不 稳定 、 求 解 精度 不 高 的 缺点 。 


改进 之 站 


江 


部 新 解 的 产生 实际 上 仅 靠 响 度 项 


1 
在 当前 全 局 最 优 粒子 周转 


整 产 生 ， 


er 
| 


有 很 大 的 随机 性 , 往往 


使 新 解 过 


因此 这 里 设计 扰动 控制 因子 a 来 


大 的 偏离 当前 全 局 最 优 解 , 或 者 调整 不 足 , 导 致 种 群 
多 样 性 的 丧失 ,从 而 产生 求解 不 稳定 和 ; 


青 度 不 高 的 缺点 。 


空 制 高 斯 扰动 项 的 干扰 
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范围 ( 式 (9)), 增强 算法 稳定 性 和 增加 种 群 多 样 性 。 关 于 a 
的 取 值 , 本 文通 过 仿真 实验 来 说 明 , 当 a 取 值 为 1、0.5、0.1、 
0.05、0.01 时 ， 分 别 用 GDEDBA 算法 对 测试 函数 Fl 和 F2 
(迭代 1 000 次 ， 维 数 为 30) 分 别 独立 进行 求解 30 次 取 其 
平均 值 ， 结 果 如 表 2 所 示 。 
表 24a 不 同 取 值 下 仿真 结果 对 比 
Table 2 Simulation results of a under different values 
a=l a =0.5 a=0.1 a =0.05 a=0.01 
Fl 3.68E-24 5.08E-62 3.19E-117 7.50E-104 3.85E-50 
F2 5.19E-16 1.06E-34 4.62E-60 2.824E-52 4.95E-26 


通过 观察 可 以 发 现 ， 在 “=0.1 时 函数 的 求解 结果 最 好 ， 


所 以 ， 本 文中 a 


经 过 以 上 分 书 ， | 本 


置 的 信息 来 保证 


取 值 为 0.1， 可 以 保 ii 


扰动 项 能 


E 算 法 的 最 优 性 能 。 
效 利 用 当前 全 局 最 优 位 
新 解 中 所 有 粒子 之 间 的 信息 交流 与 学 习 ， 


方面 可 以 有 效 帮 助 粒子 避免 陷入 局 部 最 优 ， 提 高 算法 搜索 精 


度 ， 男 一 方 


面 提高 了 种 群 多样 性 ， 加 
以 有 效 提 高 算法 的 寻 


hE 
rE 性 能 。 


2.4 


a) 设 置 相 关 


自 


(个 


算法 实现 流程 
参数 并 初始 化 蝙蝠 种 群 
永 冲 频 率 万、 位 置 X;、 


速度 Vi、 


b) 计 算 初始 化 种 群 的 函数 适应 


代 终 止 条 件 ; 
Cc) 按照 


脉冲 速率 * 和 声音 


快 了 算法 收敛 速度 ， 可 


， 随 机 生成 蝙蝠 种 群 
响 度 人 4; 


! 式 (7) 计算 指数 递减 策略 


度 值 ， 


因子 w， 


判断 是 否 达 到 和 迭 


用 式 (1) (8) 


(3) 调整 蝙 册 


日 个 体 的 频率 并 更 新 速 
位 置 x“ 并 记录 下 来 ; 


度 和 位 置 ， 


找到 当前 最 优 


d) 生 成 一 个 随机 数 rand， 并 判 册 
(9) 生成 高 斯 扰动 项 , 在 当前 迭代 次 数 的 全 局 最 优 位 置 附近 


生成 新 解 Xi ; 
©) 判断 是 否 


满足 随机 数 rand< 4 且 适 


(X”), 同时 满足 就 接受 这 个 新 的 解 ， 
最 优 解 ， 并 计算 出 Xi 对 应 的 的 适应 


(5), 
数 ; 
人 判断 fvew 


X*， 


多 判断 是 否 


是 否 为 最 优 适 应 度 值 ， 


记录 此 时 的 适应 度 值 为 目 | 


度 值 fyew 9 然 后 
(6) 更 新 5 和 4 ， 本 文 将 5 和 人 设置 为 一 个 固定 的 党 


是 否 (rand>*)， 利 用 式 


应 度 值 f (Xi 


找到 当前 全 局 最 优 解 


否 满足 算法 结束 条 


若 未 满足 ， 返 回 步骤 0); 


最 优 解 和 最 优 目 标 函 数 适 应 度 值 。 


为 了 更 加 直观 地 理 


晤 算法 的 算法 实现 流程 如 图 3 所 示 。 
3 ”仿真 实验 
3.1 算法 测试 函数 


为 了 测试 改进 算法 的 寻 优 性 能 ,选取 15 个 经 
B2425], 通过 对 基本 蝙蝠 算法 (BA)01、 标 准 
P9、 基 本 布谷 鸟 算 法 〈CS ) 2 做 横向 对 


解 改进 蝙蝠 算法 的 实现 流程 ， 


改进 旷 


测试 函数 


粒子 群 算法 (PSO ) 
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算法 指数 递减 权重 蝙蝠 算法 (IWBA) PE、 动态 调整 惯性 权 
重 的 自 适 应 蝙蝠 算法 (DAWBA) B91 作 纵向 对 比 ， 分 别 进行 


寻 优 计算 并 记录 结果 。 测 试 函 数 表 3 所 示 。 
初始 化 蝙蝠 种 群 频率 局 _ 位置 X,、 速 度 V 、 脉 冲 速 率 
声音 响 度 4 
计算 初始 种 群 的 适应 度 值 ， 并 记录 全 局 最 优 适 应 度 值 f 


生成 指数 递减 策略 因子 ， 调 整 频率 并 更 新 速度 和 位 置 ， 
记录 当前 最 优 位 置 达 ” 


否 | 随机 生成 蝙蝠 ， 计 


随机 生成 rand，rand> 和 算 其 适应 度 值 | 


时 
是 


生成 高 斯 扰动 项 ， 在 当前 迭代 次 数 的 全 局 最 优 位 
置 附近 生成 新 解 X'，， 并 计算 新 解 的 适应 度 值 


， 且 新 解 的 适应 度 值 < 当前 欠 代 最 优 适 六 


时 
是 


接受 该 新 解 X' ,并 计算 其 对 应 的 适应 度 值 hw，; 更 新 万 和 A4 
hi < hn 否 
该 新 解 为 全 局 最 优 解 , jaw 为 全 局 最 优 适应 度 值 


法 结束 条 


显示 全 局 最 优 解 和 最 优 目 标 函 数 适应 度 值 
习 3 ”改进 蝙蝠 算法 的 算法 流程 
Fig.3 Flow chart of improved bat algorithm 
表 3 列 出 了 15 个 测试 函数 的 函数 名 、 对 应 的 函数 表 i 
以 及 理论 最 优 值 。 为 了 方便 下 文 描 述 ， 对 其 进行 编号 
Fl 是 非 线 性 连续 单 峰 函 数 ， 难 以 求 到 最 优 值 ， 常 用 
法 上 ; 函数 F2 是 连续 的 平滑 多 峰 函 
变量 趋 近 于 无 穷 大 时 ， 函 数 会 形成 大 量 局 部 极 值 区 
较 高 的 寻 优 难度 ; 函数 F3 由 于 其 各 维 上 的 变量 是 密切 相关 、 
互相 作用 的 ， 函 数 解释 变量 随 着 自 量变 而 改变 ， 函 数 的 分 布 
Ph 存在 大 量 局 部 极 值 因此 很 难 取得 最 优 解 ， 可 以 检测 算法 跳 
出 局 部 的 能 力 ; 函数 F4 有 1 个 全 局 极 小 值 点 {0,0)=0， 该 函 
数 在 狭长 的 全 局 极 值 点 周围 拥有 很 多 的 局 部 极 值 ， 主 要 用 以 
测试 算法 的 收敛 性 能 ， 函 数 F5 是 一 个 复杂 的 多 峰 问题 ， 


颖 呼 攻关 
tt 于 粮 洋 此 


ul 


有 大 量 的 局 部 极 值 点 ， 主 要 用 于 检验 算法 的 群体 多 样 性 ， 函 
数 F6 是 虽然 单 峰 连续 函数 ， 具 有 一 个 全 局 极 小 点 f(1,1)=0， 


但 它 却 是 病态 的 (螺旋 型 ), 它 在 全 局 极 小 点 邻近 的 狭长 区 域 
内 取 值 变化 极为 缓慢 ， 会 发 生 振 荡 ， 可 用 于 评价 算法 的 收敛 
性 能 ; 函数 F7~F10 以 及 F14， 最 小 值 均 为 非 零 常数 ， 其 中 函 
数 F7 共有 760 个 局 部 极 小 值 ，F9 有 6 个 局 部 极 小 点 ， 常 用 
来 测试 算法 的 求解 成 功率 ; 函数 Fl11~ 函 数 13 以 及 函数 15 用 
来 测试 算法 的 求解 稳定 性 和 求解 成 功率 。 
通过 上 述 分 析 可 知 ,文中 选取 多 种 不 同类 型 的 测试 函数 ， 
以 求全 方位 测试 算法 性 能 ， 其 中 不 乏 求 解难 度 极 高 的 测试 函 
数 ， 这 里 为 了 增加 求解 难度 ， 把 所 有 函数 的 求解 维度 都 设 为 
30 维 。 
3.2 ”实验 环境 及 参数 设置 

在 仿真 实验 中 ， 仿 真 环 境 的 硬件 配置 为 Windows7 旗舰 
版 操作 系统 ，CPU 为 Intelt Core i5 2.40 GHz 处 理 器 ，6 GB 
内 存 ， 仿 真 环 境 的 软件 配置 为 MATLABR2014a。 

为 保证 比较 的 公平 性 , 这 里 对 各 个 算法 的 参数 作出 说 明 。 


比 ,， 通过 对 改进 蝙 晤 


对 于 BA、IWBA、DAWBA、GDEDBA 的 参数 设置 都 设置 为 
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图 
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ch 
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种 群 规模 n=40, 响 度 A=0.25, 脉冲 r=0.5, 最 大 频率 f=2， 
=0; 对 于 PSO， 参 数 设置 为 种 群 规模 n=40， 学 


最 小 频率 Ja 


等 : 基于 高 斯 扰动 和 指数 递减 策略 的 改进 蝙蝠 算法 


为 种 群 规模 n=40， 概 率 pa=0.25; 
高 求解 难度 ， 订 


设 为 1 000 次 ; 为 了 提 
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以 上 所 有 算法 和 狗 代 次 数 都 


E 明 算法 的 优越 性 ，15 


习 因 子 cl1=2，c2=2， 人 惯性 权重 w=0.5; 对 于 CS， 参 数 设 置 。 个 测试 函数 的 求解 维 数 都 设 为 30 维 。 
表 3 测试 函数 
Table 3 Test function 
编号 函数 名 函数 表达 式 最 优 值 
Fl Sphere OO = > 生 0 
i=l 
F2 Schwegel's Problem 2.22 £0= > x+ 工 [al 0 
i=l i=l 
1 
F3 Griewank b(n) er le -)+1 0 
F4 Ackley's f(x)=20+e—20exp(-0.2 FE Sx) cpt Sicos(27x)) 0 
na na 
FS Rastrigin's f(x)= EE —10cos(2xx,)+10] 0 
i=l 
n=-l 
F6 Rosenbrock's 万 CO = > L100 — x?) +(—x7) 0 
i=l 
F7 Hansen f= DicosG -Dn ti jcos(j+D) :tn +) -176.54 
i=1 j= 
5.1 ， 1 
F8 Branin fx) =(% 好 一 6)2+10.( )cosx +10 0.39789 
4x 8 
-6 
F9 Six-Hump Camel-Back 万 CDO =4x7 —2.1xt + 区 十 交 司 一 4 好 T4 寺 -1.0316 
folD)=+0 + + D9-140 +30 -14x +6nx +3x3)] 
ey GoldstemPnce x[30+(Cx -322)2 (18—32x +12x7 + 482 -36x 十 27 妆 )] 
F11 Tabel f=105% + 0 
i=2 
F12 Matyas fis =0.26(x? +y2) 一 0.48xy 0 
F13 Eggcrate 用 = 了 妇 +T7+T25(GSin2x+sin2 y) 0 
F14 Easom fua(X¥) = 一 cos(C)cos(OPp)exp[ 一 (26 — NA) 十 Co — Xx))] -1 
d i 
Fl15 Axis parallel hyper-ellipsoid fis(x*) = 2 0 
il j= 
3.3 寻 优 性 能 分 析 3.4 收敛 性 能 分 析 
为 了 减少 随机 因素 对 算法 求解 结果 的 影响 ， 把 每 个 函数 标 ， 算 法 的 收 


分 别 独立 运行 50 
次 的 最 优 解 (Best)、 平 均值 (Avg) 和 方差 
解 精度 最 高 的 解 进行 加 粗 ， 结 果 对 比如 表 4 
表 4 可 以 看 出 ， 
可 以 对 其 中 的 11 个 求 到 最 


优 解 ， 


次 ， 分别 计 算 15 个 测试 函数 独立 运行 50 


对 于 15 个 测试 函数 而 言 ，GDEDBA 
其 余 4 个 的 寻 优 结果 远 远 优 


于 其 他 算法 ， 说 明 改 进 算法 的 求解 达 优 率 很 高 。 仔 纪 


对 于 函数 F7~F10 和 F14, 本 文选 取 的 五 个 算法 以 及 GDEDBA 


分 析 ， 对 数 。 


都 可 以 求 到 最 优 解 ， 说 明 本 文选 取 的 所 
F7~F10 和 F14 都 显示 了 优良 的 求解 性 能 ， 


了 算法 在 求解 函数 
但 是 对 于 函数 F7 
和 Fl10 来 说 , GDEDBA 的 求解 结果 方差 最 小 , 说明 改进 后 的 曲线 几 


速 下 降 , 万 


算法 的 收敛 速度 是 衡量 算法 性 能 的 重要 指 
皇 


敛 图 可 以 更 加 直观 地 对 比 各 种 算法 的 收敛 性 
(Var)， 并 把 求 示 了 BA、PSO、CS、 
所 示 。 优化 算法 在 求解 15 个 标 ;# 

图 中 横 轴 


能 。 图 4~18 展 


IWBA、DAWBA 和 GDEDBA 这 6 种 
侍 测 试 函 数 时 的 适应 度 值 收敛 曲线 。 
为 迭代 次 数 ， 纵 轴 为 适应 度 值 。 男 外 ， 为 了 更 清楚 
地 显示 收敛 情况 ， 把 函数 FlI~F6，F11~F13 和 F15 的 纵 轴 取 


图 4~8、 图 14~ 16 以 及 图 18 可 知 ,在 求解 函数 Fl1~F5、 
Fl11~F13 以 及 F15 时 , GDEDBA 的 收敛 


线 向 图 的 右 下 角 迅 


其 是 图 6~8 对 应 的 函数 F3~FS，GDEDBA 的 收敛 
乎 接近 垂直 下 降 ， 收 敛 非常 快 ， 


于 其 余 算 法 ; 


远 远 高 


算法 稳定 性 更 强 ; 对 于 函数 F4~F6, PSO, IWBA 和 GDEDBA ” 观察 图 6~8 可 发 现 ， 其 余 5 种 算法 在 收敛 过 程 中 出 现 了 不 同 
可 以 求 到 最 优 解 ， 而 其 余 算法 求 不 到 最 优 解 ， 对 于 函数 程度 的 阶梯 状 , 观 察 图 14~16 发 现 ,PSO 在 求解 函数 F11~F13 
F1,F2,F11 和 F15 本 文选 取 的 五 个 算法 以 及 改进 后 的 算法 都 时， 在 欠 代 前 期 出 现 了 剧烈 震荡 ， 而 GDEDBA 收敛 曲线 平 
没有 求 到 最 优 解 ， 但 是 GDEDBA 的 求解 精度 高 出 其 余 算法 滑 ， 由 收敛 曲线 可 以 看 出 ，GDEDBA 求解 稳定 性 很 强 ， 显示 
近 百 倍 的 数量 级 ， 而 且 求 解 方 差 远 远 高 于 其 余 算 法 ， 说 明 了 很 强 的 鲁 棒 性 。 由 图 10~13 以 及 图 17 可 知 ， 求 解 函数 
GDEDBA 的 稳定 性 更 强 ， 而 对 于 函数 F3,F12 和 F13， 只 有 F7~F10 以 及 F14 时 ，GDEDBA 的 收敛 曲线 几乎 与 其 余 算法 
本 文 的 改进 算法 GDEDBA 可 以 求 到 最 优 解 ， 其 余 算法 都 没 ” 重合 ， 上 文 可 知 F7~F10 以 及 F14 都 可 以 求 到 最 优 解 ， 

有 求 到 最 优 解 。 经 过 以 上 分 析 可 以 得 出 ， 本 文 设计 的 改进 算 两 者 互 为 证 明 ， 算 法 的 收 钱 速度 与 收敛 精度 是 存在 一 定 联系 
法 在 在 寻 优 性 能 上 明显 优 于 其 他 算法 ， 显 示 了 更 高 的 寻 优 精 的 ,由 图 9 可 知 ,PSO 的 收敛 速度 快 于 其 余 算法 ,而 DGEDBA 


度 和 更 高 的 稳定 性 。 


和 IWBA 的 收敛 曲线 呈现 重合 趋势 ，] 
数 F6，PSO、IWBA 和 GDEDBA 都 可 


由 上 文 得 知 ， 对 于 函 
以 求 到 最 优 解 。 
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录用 定稿 
表 4 ， 寻 优 精度 的 对 比 结果 
Table 4 Comparison results of optimization accuracy 
函数 算法 Best Avg Var 
BA 9.0151E-04 1.4130E-03 3.4986E-08 
PSO 2.0146E-01 5.3673E-01 3.5194E-02 
CS 1.6948E-05 5.1693E-05 7.2214E-10 
IWBA 6.0092E-05 2.8104E-04 8.0002E-09 
DAWBA 3.7525E-07 1.1544E-06 3.4034E-13 
GDEDBA 4.4294E-123 2.1932E-117 5.1902E-233 
BA 1.2693E-01 2.1723E-01 1.78E-02 
PSO 1.309 2.680037255 0.425097685 
CS 1.2434E-01 2.6981E-01 6.1909E-03 
IWBA 7.5163E-02 3.2902E-01 6.10E-02 
DAWBA 3.4754E-03 1.3185E-01 3.92E-02 
GDEDBA 3.9700E-63 3.426SE-60 S.88S1E-119 
BA 3.6417E-05 7.2499E-05 .13885E-10 
PSO 1.0912E-02 2.2616E-02 6.3948E-05 
CS 5.1171E-05 7.0376E-04 1.2800E-06 
IWBA 8.7483E-06 3.0170E-05 .19875E-10 
DAWBA 1.9098E-08 8.2571E-08 2.42595E-15 
GDEDBA 0 0 0 
BA 4.9913E-05 2.9146E-04 .64453E-08 
PSO 0 0 0 
CS 0 0 0 
F4 
IWBA 0 0 0 
DAWBA 9.7607E-08 2.0547E-06 1.96856E-12 
GDEDBA 0 0 0 
BA 1.78E-15 1.31E-10 4.97312E-20 
PSO 0 0 0 
CS 0 0 0 
FS 
IWBA 0 0 0 
DAWBA 0 6.77E-14 1.7116E-26 
GDEDBA 0 0 0 
BA 1.49E-09 8.96E-08 8.48199E-15 
PSO 0 0 0 
CS 0 1.1434E-25 1.4415E-49 
Ee IWBA 0 0 0 
DAWBA 5.25E-13 9.11E-11 2.54143E-20 
GDEDBA 0 0 0 
BA -176.54 -171.1060784 176.2489323 
PSO -176.54 -172.9331373 135.213922 
CS -176.54 -176.36676 85.762971 
Ey IWBA -176.54 -171.7623529 196.5068224 
DAWBA -176.54 -174.7129412 54.47850918 
GDEDBA -176.54 -175.3219608 37.07565208 
BA 0.39789 0.39789 5.02899E-32 
PSO 0.39789 0.39789 5.02899E-32 
CS 0.39789 0.39789 2.87239E-32 
四 IWBA 0.39789 0.39789 5.02899E-32 
DAWBA 0.39789 0.39789 5.02899E-32 
GDEDBA 0.39789 0.39789 5.02899E-32 
BA -1.0316 -1.0316 4.52609E-31 
PSO -1.0316 -1.0316 4.52609E-31 
CS -1.0316 -1.0316 4.59582E-31 
本 IWBA -1.0316 -1.0316 4.52609E-31 
DAWBA -1.0316 -1.0316 4.52609E-31 
GDEDBA -1.0316 -1.0316 4.52609E-31 
BA 4.058823529 28.01647059 
PSO 3 0 
CS 3 0 
F10 
IWBA 3 4.588235294 41.16705882 
DAWBA 3 4.058823529 28.01647059 
GDEDBA 3 3 0 


函数 算法 Best Avg Var 
BA 3.5727E-04 6.9015E-01 0.477236782 
PSO 1.8816E-123 2.6125E-110 1.9336E-218 
| CS 2.1420E-25 9.4176E-22 5.7827E-42 
IWBA 1.1134E-51 7.9783E-46 5.51478E-90 
DAWBA 3.5727E-04 6.9015E-01 0.477236782 
GDEDBA 1.9677E-270 6.8276E-258 0 
BA 1.3617E-11 8.8821E-10 6.2287E-19 
PSO 7.0720E-161 3.1377E-148 2.1632E-294 
i CS 1.3575E-45 1.6298E-40 5.3858E-79 
IWBA 4.7641E-76 2.9422E-69 3.1479E-136 
DAWBA 1.3617E-11 8.8821E-10 6.2287E-19 
GDEDBA 0 0 0 
BA 2.7643E-09 2.3366E-07 5.09804E-14 
PSO 2.3539E-158 1.5851E-147 1.2409E-292 
二 CS 1.1081E-39 1.7153E-34 5.3844E-67 
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F14 
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i CS 2.0909E-04 5.9789E-04 8.0044E-08 
IWBA 5.5066E-03 2.1726E-02 1.5174E-04 
DAWBA 6.7472E-06 2.3829E-05 1.3612E-10 
GDEDBA 1.1240E-121 6.3033E-117 2.7792E-232 
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Fig.5 Convergence curves of function F2 
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图 12 函数 F9 的 收敛 曲线 
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图 15 函数 F12 的 收敛 曲线 


Fig. 15 Convergence curves of function F12 
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函数 F15 的 收敛 
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图 13 函数 F10 的 收敛 曲线 


Fig. 13 Convergence curves of function F10 
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图 14 函数 Fll 的 收敛 曲线 


Fig. 14 Convergence curves of function Fl11 
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图 17 函数 Fl4 的 收敛 曲线 


Fig. 17 Convergence curves of function F14 


综合 以 上 分 析 ,， GDEDBA 在 处 理 单 峰 多 维 、 多 峰 多 维 的 


收敛 速度 和 求解 稳定 性 方面 比 其 余 五 种 算法 


求解 性 能 优良 。 
4 ”结束 语 


针对 基本 蝙蝠 算法 后 中 


问题 上 可 以 显示 很 高 的 收敛 速度 与 收敛 精度 ， 并 且 稳 定性 
于 其 他 算法 , 表明 了 改进 后 的 算法 (GDEDBA ) 在 寻 优 精度 、 
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有 显著 的 优越 性 ， 


收敛 速度 慢 、 寻 优 精度 低 、 求 解 
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Fig. 18 ”Convergence curves of function F15 


不 稳定 等 不 足 , 通过 在 产生 局 部 新 角 
方差 为 当前 和 迭代 次 数 的 最 优 解 绝对 
长 控制 因子 确保 扰动 幅度 ， 既 提高 


Ba 


坚 的 公式 中 加 入 均值 为 0、 
值 的 高 斯 扰动 项 ， 加 入 步 
了 算法 的 精确 搜索 能 力 ， 
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又 增加 种 群 多 样 性 与 活跃 性 ， 进 一 步 挖掘 蝙蝠 算法 的 仿生 性 
能 ， 对 速度 更 新 公式 加 入 指数 递减 策略 使 算法 迅速 进入 局 部 
精确 搜索 ， 实 现 对 基本 蝙蝠 算法 从 全 局 寻 优 到 局 部 寻 优 都 做 
了 有 效 改进 来 提高 算法 性 能 。 仿真 结果 表明 , GDEDBA 的 寻 
优 精度 ,收敛 速度 以 及 算法 稳定 性 方面 均 有 很 大 程度 的 提升 ， 
寻 优 性 能 优良 。 


参考 文献 : 


[1] Yang Xinshe. A new meta-heuristic bat-inspired algorithm [Cl]// Proc of 


Nature Inspired Cooperative Strategies for Optimization. Berlin 
Eidelberg: Springer-Verlag, 2010: 65-74. 

[2] TsaiP W, Pan J S, Liao B Y, et al. Bat algorithm inspired algorithm for 
solving numerical optimization problems [J]. Applied Mechanics and 
Materials, 2012, 148-149 (148): 134-137. 

[3] 李煜, 装 宇 航 , 刘 景 厅 , 融合 均匀 变异 与 高 斯 变异 的 蝙蝠 优化 算法 
四 . 控制 与 决策 , 2017, 32 (10): 1775-1781. (Li Yu, Pei Yuhang, Liu 
Jingsen. Bat optimal algorithm combined uniform mutation with 
Gaussian mutation [J]. Control and Decision, 2017, 32 (10): 
1775-1781.) 

[4] Yang Xinshe. Bat algorithm for multi-objective optimization [J]. 
International Journal of Bio-Inspired Computation, 2011, 3 (5): 
267-274. 

[5] Kang M, Kim J, Kim J M, et al. Reliable fault diagnosis for incipient 
low-speed bearings using fault feature analysis based on a binary bat 
algorithm [J]. Information Sciences, 2015, 294 (2): 423-438. 

[6] Abd-Elazim S M, Ali E S. Load frequency controller design via BAT 
algorithm for nonlinear interconnected power system [J]. Electrical 
Power and Energy Systems, 2016, 77 (5): 166-177. 

[7] Singh K, Vasant P. Elamvazuthi I, et al. PID tuning of servo motor 
using bat algorithm [J]. Procedia Computer Science, 2015, 60 (1): 
1798-1808. 

[8] 马 邦 雄 , 叶 春 明 . 基于 蝙蝠 退火 算法 的 无 等 待 流水 线 调度 问题 研究 
[]. 数学 理论 与 应 用 ，2014，34 (1): 92-101. (Ma Bangxiong，Ye 
Chunming. A bat annealing algorithm for solving the no-wait flow shop 
scheduling problem [J]. Mathematical Theory and Applications, 2014, 
34 (1): 92-101.) 

[9] 金 伟 健 , 王 春 枝 . 基于 蝙蝠 算法 的 云 计 算 资 源 分 配 研究 [J]. 计算 机 
应 用 研究 , 2015, 32 (4): 1184-1187. (Jin Weijian, Wang Chunzhi. Study 
on bat algorithm in cloud computing resources allocation [J]]. 
Application Research of Computers, 2015, 32 (4): 1184-1187. ) 

[10] Ye Zhiwei, Wang Mingwei, Liu Wei, et al. Fuzzy entropy based optimal 
thresholding using bat algorithm [J]. Applied Soft Computing, 2015, 31 
(6): 381-395. 

[11] 陈 海 手 . 改进 蝙蝠 算法 优化 极限 学 习 机 的 图 像 分 类 [加 . 激光 杂志 ， 
2014, 35 (11): 26-29. (Chen Haiting. Image classification based on 
extreme learning machine optimized by improved bat algorithm [J]. 
Laser Journal, 2014, 35 (11): 26-29. ) 

[12] Selim Y, Ecir U. Kiiciiksille. A new modification approach on bat 
algorithm for solving optimization problems [J]. Applied Soft 
Computing, 2015, 28 (12): 259-275. 

[13] Meng Xianbing, Liu Yu, Zhang Hengzhen, et al. A novel bat algorithm 
with habitat selection and doppler effect in echoes for optimization [J]. 
Expert Systems with Applications, 2015, 42 (17): 6350-6364. 

[14] Luo Jun, Liu Liheng, Wu Xianyi. A double-subpopulation variant of the 


bat algorithm [J]. Applied Mathematics and Computation, 2015, 263 (7): 


宋 一 民 ， 等 : 基于 高 斯 扰动 和 指数 递减 策略 的 改进 蝙蝠 算法 


361-377. 

[15] Nakamura R Y M, Pereira L, Costa K, et al. BBA: a binary bat 
algorithm for feature selection [C]// Proc of the 25th SIBGRAPI 
Conference on Graphics, Patterns and Images. Piscataway: IEEE Press, 
2012: 291-297. 

[16] He Xingshi, Ding Wenjing, Yang Xinshe. Bat algorithm based on 
simulated annealing and Gaussian perturbations [J]. Neural Computing 
and Applications, 2013, 25 (2): 459-468. 

[17] 成 远 航 , 蔡 延 光 , 蔡 侨 ,等 . 旅行 商 问 题 的 混沌 混合 离散 蝙 蝙 算法 
[四 . 电子 学 报 , 2016, 44 (10): 2543-2547. (Qi Yuanhang, Cai Yanguang， 
Cai Hao, et al. Chaotic hybrid discrete bat algorithm for traveling 
salesman problem [J]. Acta Electronica Sinica, 2016, 44 (10): 
2543-2547. ) 

[18] Premkumar K, Manikandan B V, Speed control of Brushless DC motor 
using bat algorithm optimized adaptive neuro-fuzzy inference system 
[J]. Applied Soft Computing, 2015, 32 (7): 403-419. 

[19] Shi Yuhui, Eberhart R C. Empirical study of particle swarm 
optimization [Cl]// Proc of Congress on Evolutionary Computation. 
Washington DC: IEEE Press, 1999: 1945-1949. 

[20] Shi Yuhui, Eberhart R C. Fuzzy adaptive particle Swarm optimization 
[CI]/ Proc of Congress on Evolutionary Computation. Piscataway,NJ: 
IEEE Press, 2001: 101-106. 

[21] Eberhart R C, Shi Yuhui. Tracking and optimizing dynamic systems 
with particle swarms [Cl]// Proc of Congress on Evolutionary 
Computation. Piscataway,NJ: IEEE Press, 2001: 101-106. 

[22] 陈 贵 敏 , 机 建 捷 ， 韩 琪 . 粒子 群 优化 算法 的 惯性 权 值 递减 策略 研究 
四 . 西安 交通 大 学 党报,，2006, 40 (1): 53-61. (Chen Guimin, Jia 
Jianyuan, Han Qi. Study on the strategy of decreasing inertia weight in 
particle swarm optimization algorithm [J]. Journal of Xi”an Jiaotong 
University, 2006, 40 (1): 53-61. ) 

[23] 朱德 刚 ， 孙 涟 ,起 嘉 ， 等 . 基于 高 斯 扰动 的 粒子 人 群 优化 算法 [加 计 
算 机 应 用 , 2014, 34 (3): 754-759. (Zhu Degang, Sun Hui, Zhao Jia, et 
al. Particle swarm optimization algorithm based on Gaussian 
disturbance [J]. Journal of Computer Applications, 2014, 34 (3): 
754-759. ) 

[24] Jordehi A R. Chaotic bat swarm optimization (CBSO) [J]. Applied Soft 
Computing, 2014, 26 (1): 523—530. 

[25] 尚 俊 娜 ， 刘 春 菊 , 后 克 强 , 等 具有 自学 习 能 力 的 变异 蝙蝠 优化 算法 
及 性 能 仿真 [机 . 系统 仿真 学 报 , 2017, 29 (2): 301-307. (Shang Junna， 
Liu Chunju, Yue Keqiang, et al. Variation bat algorithm with 
self-learning capability and its property analysis [J]. Journal of System 
Simulation, 2017, 29 (2): 301-307. ) 

[26] Kennedy J, Eberhart R. Particle swarm optimization [C]/ Proc of IEEE 
International Conference on Neural Networks. Piscataway,NJ: IEEE 
Press, 1995: 1942-1948. 

[27] Yang Xinshe, Deb S. Cuckoo search via Lévy flights [Cl]// Proc of 
World Congress on Nature & Biologically Inspired Computing. 
Piscataway,N]J: IEEE Press, 2009: 210-214. 

[28] 宋 一 民 . 基于 改进 蝙蝠 算法 的 物流 中 心 选 址 问题 研究 [D]. 开 封 : 河 
南大 学 , 2017. (Song Yimin. Research on logistics center location based 
on improved bat algorithm [D]. Kaifeng: Henan University, 2017. ) 

[29] 甘 字 航 , 刘 景 杰 ,李煜 . 一 种 动态 调整 惯性 权重 的 自 适应 蝙蝠 算法 
四 . 计算 机 科学 , 2017, 44 (6): 240-244. (Pei Yuhang, Liu Jingsen, Li 
Yu. An adaptive bat algorithm with dynamically adjusting inertia weight 
[J]. Computer Science, 2017, 44 (6): 240-244. ) 


